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Abstract. The question of how the scattering cross section changes when the spectra of 
the colliding nuclei have low-excitation particle-emitting resonances is explored using a multi- 
channel algebraic scattering (MCAS) method. As a test case, the light-mass nuclear target 8 Be, 
being particle-unstable, has been considered. Nucleon-nucleus scattering cross sections, as well 
as the spectra of the compound nuclei formed, have been determined from calculations that 
do, and do not, consider particle emission widths of the target nuclear states. The resonant 
character of the unstable excited states introduces a problem because the low-energy tails 
of these resonances can intrude into the sub-threshold, bound-state region. This unphysical 
behaviour needs to be corrected by modifying, in an energy-dependent way, the shape of the 
target resonances from the usual Lorentzian one. The resonance function must smoothly reach 
zero at the elastic threshold. Ways of achieving this condition are explored in this paper. 



1. Introduction 

The advent of radioactive ion beam (RIB) physics prompts the consideration of new theoretical 
challenges involving weakly-bound systems. Radioactive nuclei, especially those close to the 
drip lines, can have quite low particle emission thresholds and consequently have low- lying 
resonance states in their spectra. Low-energy scattering of RIB is a coupled-channel problem 
that involves such low-lying resonant states of the scattered ion. To address this, a Multi- 
Channel Algebraic Scattering (MCAS) formalism [I] is used, in which momentum space solutions 
of coupled Lippmann-Schwinger equations are found, including also all negative-energy bound 
states. A finite-rank separable representation using an "optimal" set of sturmian functions [2J 
serves to construct an input matrix of nucleon-nucleus interactions. The MCAS method has the 
ability to locate all compound system resonance centroids and widths, regardless of how narrow. 
For full details see reference [I]. Also, by use of orthogonalizing pseudo-potentials (OPP) in 
generating sturmians, it ensures the Pauli principle is not violated [3], despite the collective 



model formulation of nucleon-nucleus interactions used. Otherwise, some compound nucleus 
wave functions possess spurious components. 

In the following sections, we first briefly summarize the MCAS method [HE], which has been 
extended to unstable target states in Ref. [4j. In section II, we display and discuss our results 
for neutron scattering from 8 Be, with and without the excited target states having non-zero 
widths. In section III, we look briefly at the implications of letting the target-state widths have 
energy dependence, to ensure that these resonances do not extend to negative energies, where 
they would produce unphysical behaviour for bound, sub-threshold, states. Finally, section IV 
gives our brief concluding remarks. 

2. The MCAS formulation for unstable target states 

S'-matrix equations in the MCAS methodology take the form: 
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Traditionally, all target states are taken to have eigenvalues of zero width. Then the integrals 
in the (complex) Green's functions are evaluated using the method of principal parts, where, in 
the limit e — > 0, the Green's functions take the form 
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rj is a column vector of sturmian eigenvalues and x are form factors determined from the chosen 
sturmian functions. The factor \i in Eqs. [J and [2] is \i = 2m re( i/h 2 , with m re d being the reduced 
mass. This method assumes time evolution of target states is given by 

|a;, t) = e -' lHot,n \x, t ) = e - iEot/h \x, t ) ■ (3) 

However, if states decay, they evolve as 

^ t)=e -(Tt/2H) e -iE t/H (4) 

Thus, in the Green's function, channel energies become complex, as do the squared channel wave 
numbers, 
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where ?f is half the width of the target state associated with channel c. The Green's function 
matrix elements are then 
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When poles are moved significantly off the real momentum axis, integration along this axis is 
feasible. 



3. The case of 9 Be 

The low excitation spectrum of 8 Be has a + ground state that decays into two a-particles 
(width: 6 x 1(T 6 MeV), a 2+ resonance state (centroid: 3.03 MeV; width: 1.5 MeV) and a 4+ 
resonance state (centroid: 11.35 MeV; width: ~3.5 MeV). In Table [Tj are given the parameters 



Table 1. The parameter values defining the nucleon- 8 Be intera ction 





Odd parity 


Even parity 


Central (MeV) 


-31.5 


-42.2 


V u (MeV) 


2.0 


0.0 


Vis (MeV) 


12.0 


11.0 


V„ (MeV) 


-2.0 


0.0 


Geometry Rq 


= 2.7 fm a = 


0.65 fm fa = 0.7 




Osi/2 


0p 3 /2 


o+ \(opn 


1000 


0.0 


2+ A(° pp ) 


2.0 


0.2 


4 + x (opp) 


0.0 


0.0 



in the channel-coupling potentials, as well as the A values giving the strengths of Pauli 
hindrance. Two evaluations of the n+ 8 Be cross section are obtained using this spectrum; one 
with target-state widths set to zero, and the other taking into account the widths of the excited 
levels. The two results will be identified by the terms 'no- width' and 'width', respectively. Both 
cases are calculated with the same nuclear interaction, taken from a rotor model [5]. Tabled 
displays 9 Be spectrum results from these calculations. 

Table 2. Widths of resonances in n+ 8 Be scattering for the no-width and width cases, compared 
to experimental values. Boldface entries highlight the matches we find between theory and 
experiment to within a factor of 3. 





r 

1 exp. 


^no— width 


* no— width 


r width 


^ width 
r exp. 


1+ 

2 


0.217±0.001 






1 


595 


7.350 


7~ 
2 


1.210±0.230 


2.08xl0~ 5 


1.72xl0~ 5 


1 


641 
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1.080±0.110 


0.495 


0.458 


1 


686 


1.561 


5 + 
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0.282±0.011 


0.187 


0.663 





740 


2.624 


3~ 
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1.330±0.360 


0.466 


0.350 


3 


109 


2.337 


5~ 
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7.8xl0~ 4 


0.060 


76.74 


2 


772 


3554 


9 + 
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1.330±0.090 


0.386 


0.290 


2 


498 


1.878 


3 + 
2 


0.743±0.055 


3.286 


4.423 


5 


162 


6.947 



Taking the excited states of Be to be resonances gives the same spectral list as when they 
are treated as zero- width, but the evaluated widths of the compound nuclear states significantly 
increase, as reflected in the cross sections. These increases bring the theoretical 9 Be state 
widths closer, often significantly, to experimental values, as evidenced by comparison of the final 
columns of the two tables. Only the | , also poorly recreated in centroid, has a worse match 
with widths applied. The resultant elastic and reaction cross sections are shown in figure [TJ The 
blue line depicts the results without widths, the red line with constant widths multiplied by the 
Heaviside (step) function: H(E) = for E < 0.0; H(E) = 1 for E > 0. (For the green line, see 
section HI) 

Introducing target state widths, the resonances in 9 Be are suppressed but still present; their 




Figure 1. Elastic (left panel) and reaction (right panel) cross section as a function of neutron 
laboratory energy with three treatments of the widths of target states (see text). 



widths increasing and magnitudes decreasing so as all but the , and arguably (|) 2 , cannot 
be discerned from the background. The reaction cross section found in the no-width case is zero 
until 3.4 MeV (3.03 MeV in the centre of mass frame), the first inelastic threshold. When the 
target state has width different from zero the elastic cross section has flux loss from zero projectile 
energy upwards. There is an asymptotic behaviour as energy approaches zero. In the reaction 
cross section from the width calculation, there is a broad peak at 1.5 MeV, corresponding to the 
| + resonance. Using widths not varying with energy (see the next section), the bound states, 
E < 0.0, of the compound 9 Be system are unstable due to the couplings to the decaying excited 
states of the 8 Be core. For this reason the widths have been multiplied by a Heaviside function 
to ensure that the compound system is stable at negative energies. An improvement of this 
technique is discussed in the next section. 



4. Energy dependence of widths 

A smooth description of the energy dependence of the decay width is needed to avoid any 
pathological behaviour caused by the Heaviside function. This function allows maintaining the 
stability of the compound system in the sub-threshold region, by setting to zero all widths for 
E < 0, but can generate "ghost" levels where the application of widths could move a sub- 
threshold state above E = 0, thus creating two copies of the same state. Also, a dramatic 
overestimation of the reaction cross-section occurs in the region close to threshold. 

We therefore need to improve the model with a suitable parametrization of the energy- 
dependence of the decaying width for the excited levels, with the constraint that the width 
must be zero at the scattering threshold, and equal to the experimental width at the resonance 
centroid energy. We are presently exploring various forms of energy dependences for the widths 
of these excited states (T C (E) = T c x U{E)). Examples of such expressions are 

U W = (1 + I/-1 { A {i-xy + i - A(x*) + 1 } H[E) (7) 



where x = E/E r , H(E) is the Heaviside function defined above and A is a parameter, and the 
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In this paper we study only the first form. The second, and others, will be explored in future 
work. Also yet to be taken into account is causality, as needed when dealing with energy- 
dependent widths [6]. 

In figure[2]we give the induced energy dependence for the function U(x = E/E r ) with respect 
to the renormalized energy x. The curves differ for different values of the parameter A. The 
green line in figure [T] shows the effect of this energy-dependent width on elastic and reaction 
cross sections, respectively, when A = 0.5. Clearly, for this value of A, the asymptotic behaviour 
in the reaction cross section at E = 0.0 is suppressed, but not resolved. The use of an energy 
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Figure 2. The scaling function, U(E/E r ), for various values of the parameter A. 

dependence in the widths of the excited states of the target introduces a new aspect that has 
to be taken into consideration. In principle, an energy dependence in the widths induces a 
corresponding energy-dependence in the positioning of the resonance centroids. This additional 
effect is a direct consequence of the principle of causality which demands that the real and 
imaginary part of any Green's function matrix elements are related by a dispersion integral. 
Such an aspect has been extensively discussed in the study of the properties of the nuclear 
optical potential [6]. Eq. [7] is particularly convenient in this regard, because the corresponding 
dispersion integral leads to an analytical energy-dependent shift of the resonance centroids. Such 
analytical forms have been derived in Ref. [7] for a class of functions which includes Eq. [JJ The 
resulting energy shift function AS C (E) = T c x S(E) is given by the analytical expression: 
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(1 + A) 2 m V 1 + ^ 
[l + A) 2 -l vr 



C2= (1 + ^)2 -l—{2 +arCt9VA 



(10) 



and 
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5. Conclusions 

An extension to the multi-channel algebraic scattering (MCAS) formalism that considers 
particle-decay widths of target nucleus eigenstates has been applied to a range of light mass 
nuclear targets [5], in addition to the results for n+ 8 Be shown in this paper. To ensure that any 
resonance aspect of target states does not have influence below the nucleon-nucleus threshold in 
our formalism, we have used a Heaviside function to cut off tails of the positive-energy resonances 
below zero energy of the compound system. This procedure, however, is too simplistic since it 
introduces non-physical singularities in cross sections at the approach to zero energy, and may 
also incorrectly overestimate reaction cross just above the threshold. 

A better energy dependent scaling factor is needed: one with value one at the centroid 
energy and approaching zero for lower and higher energies. We propose here possible means 
of achieving such modified Lorentzians, and show, for one example, their effect on the elastic 
and reaction cross sections in the scattering of neutrons from 8 Be. Work is in progress to apply 
these and other forms to the light-nuclear systems considered in our previous work. Work is in 
progress, also, for preserving causality in the presence of energy-dependent widths by the use of 
appropriate dispersion relations. 
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